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ABSTRACT 

We use a gravitational bar torque method to compare bar strengths (the maximum 
tangential force normalized by radial force) in B and i/-band images of 152 galaxies 
from the Ohio State University Bright Spiral Galaxy Survey. Our main motivation is to 
check how much the difference in the rest-frame wavelength could affect comparisons 
of bar strengths in low and high redshift observations. Between these two bands we 
find an average bar strength ratio Qb/h= 1-25 which factor is nearly independent 
of the morphological type. We show that Qb/h> 1 is mostly due to reduced bulge 
dilution of radial forces in the B-band. The bar torque method needs an estimate for 
the vertical scale height of the galaxy, based on the radial scale length of the disk 
and the galaxy's morphological type. Since these two might not always be possible to 
determine at high redshifts in a reliable manner, we also checked that similar results 
are obtained with vertical scale heights estimated from the radii corresponding to the 
K-hand surface brightness of 20 mag/arcsec^. Also, we made a simple test of the 
usability of the bar torque method at high redshifts by checking the effects of image 
degradation (nearest neighbour sampling without any adjustment of noise levels): we 
found that the estimated bar strengths varied by ±10% at most as long as the total 
extent of the bar was at least 10 pixels. Overall, we show that the gravitational bar 
torque method should provide a proficient tool for quantifying bar strengths also at 
high redshifts. 
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1 INTRODUCTION 

Bars are important drivers of galaxy evolution, which high- 
lights the importance of studying bar fractions and the 
properties of bars at a large range of redshifts. Bars gen- 
erally have fairly old stellar populations, thus being well 
visible in the near-IR, where also the obscuration by dust 
is minimal. However, the near-IR images we obtain from 
the distant galaxies are redshifted so that at redshift z ~ 1 
their rest-frame wavelength lies in the optical. Therefore, 
it is important to compare the optical and near-IR prop- 
erties of bars in local galaxies. For that purpose the B 
and //-band images are very convenient, offering a wide 
separation in wave lengths. Namely, bars are still promi- 
nent in the 5-band psk ridge et al.l (2002l : IWhvte et aLll2002l : 
iMenendez-Delmestre et al. 2007), where they can be iden- 
tified in a fairly reliable manner, whereas, due to the 
Balmer break, ba rs typically disappear in the ultraviolet 
ISheth et al.ll2003l 'l. On the other hand, the //-band is ah 
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ready fairly free of the effects of dust and thus offers a reli- 
able tracer for the stellar mass distribution. 

If bars were formed due to instabilities in the disk, one 
would expect a large fraction of barred galaxies at high red- 
shifts in the framework of the hierarchical clustering models. 
This is because these models ind icate dynamically colder 
disks in the distant Universe ({Navarro. Frenk fc Whit^ 
Il995h . Also, interactions of disk galaxies were probably 
more frequent in the past (jFerguson. Dickinson fc Williama 
|2000| . and reference s therein), promoting the formatio n 
of bars (discussed in lElme green. Elmegreen fc Hirstll2004l ). 
Therefore, it was surprising to find a lack of bars 
in g alaxies at redshift of z > 0.7 ( Abraham et ahl 
Ivan den Bergh et al.lll996l: Ivan den Bergh et al.ll20oH: 



1999 



van den Bergh. Cohen fc Crabbj|200ll ). Since this might be 
an artifact due to n oise, poorer image re s olution , and red- 
shifted wavelength, Ivan den Bergh et al.1 l|2002l) degraded 
about 100 /3-band images of the Ohio State University 
Bright Spiral Galaxy Survey (osuBSGS) to study these ef- 
fects. By mimicking the characteristics of the images in the 
Hubble Deep Field, using the Wide Field Planetary Camera 
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2 on HST, they created artificial images of tliose galaxies as 
they would appear in the /-band at a redshift z — 0.7. They 
found that two-thirds of bars initially classified as strong 
bars (SB), were still detectable in the degraded images, 
whereas weak bars (SAB) would largely disappear. They 
concluded that while selection effects reduce the number of 
observed bars, they cannot completely explain the lack of 
bars at high redshift. 

The ear l y stud ies of bar frequencies were challenged by 
ISheth et all (|2003l ). who found similar bar frequencies at 
high and at low redshifts for very large bars. However, the 
number statistics were a severe concern, since only four bars 
were identified. These authors used galaxies in the Hubble 
Deep Field-North, observed at V, I and //-bands, paying 
particular attention to the resolution limits due to which 
small bars canno t be detected at high redshifts. In two sub - 
sequent studies l|Elmegreen et al.l [20041 : iJogee et al]l2004 ). 
higher resolution images were used, based on the optical 
Advanced Camera for Surveys (ACS) at the Hubble Space 
Telescope. Both studies found a cons t ant ba r fraction in the 
redshift range z = - l. ljoeee et all (|2004 ) also found that 
it holds separately for the distinct intervals z = — 0.7 
and 2 = 0.7 — 1.0, using broad V and z-bands to identify 
the bars. Taking into account the wavelengths used in their 
studies this constant bar fraction is not completely expected. 
Namely, if the true bar fraction is the same at all redshifts 
and bars are more prominent in the near-IR, one would ex- 
pect a larger number of detected bars in the intermediate 
redshifts. There the rest-frame band shifting is not yet im- 
portant j_so_tlm^_both_stron^^ should be de- 
tected (|Menendez-Delmestre et al.ll2007l ). However, in the 
above studies the number statistics are still too poor (less 
than 300 galaxies spread over the entire redshift range) for 
any reliable discussion of how the bar fraction evolves as a 
function of redshift. Also, using ACS very deep images are 
required for detecting weak bars at high redshift. 

The ab o ve stu dies have been superseded recently by 
ISheth et all (|2007l ). using wide field COSMOS observations 
in /-band at z < 0.84, based on deeper images than used in 
the previous studies. Contrary to most previous studies they 
found that the fraction of bars declines rapidly with redshift. 
This drop of bars was also found to be more dramatic for 
strong bars. These studies of bar frequencies clearly show, 
besides the difficulty of measuring bar frequencies at high 
redshifts, also the importance of measuring the strengths of 
the bars at high redshifts in a reliable manner. 

In this study we use the OSUBSGS to compare bar 
strengths derived from the B and //-band images of nearby 
galaxies, using a gravitational bar torque method. Besides 
bar strength, our method gives simultaneously also the 
length and the relative luminosity of the bar. The outline 
of this paper is as follows. In Sections [2] and [3] we describe 
the details of our sample and the method used, and also ad- 
dress the factors which might affect the bar strength mea- 
surements at higher redshifts. We then show comparisons 
between B and //-band bar strengths in Section 3] the re- 
sults are discussed in Section (5] and summarized in Section 

m 



2 SAMPLE 

We use the B and //-band images from the OSUBSGS 
(|Eskridge et al.l [20o3 ). The original OSUBSGS sample con- 
sists of 205 spiral galaxies, selec ted from the Third Referenc e 
Catalogue of Bright Galaxies (|de Vaucouleurs et al.|[l99ll ). 
with the following criteria: Hubble type index ^5 T ^ 9, 
B magnitude m_g ^ 12, galaxy diameter D ^ 6'.5 and a 
declination in the range —80° ^5 (5 ^5 +50°. In this study we 
use a subsample of the OSUBSGS, with a restriction in incli- 
nation set to i < 65°, limiting the number to 156 galaxies. 
Th e properties of this sample have been discussed in detail 
bv iButa. Laurikainen fc Said l|2004h . The requirement that 
both B and //-band images must be available restricts our 
final number of galaxies to 152. Additionally, we leave out 
NGC 3338 in some of our calculations because of a bright 
saturated foreground star in the /3-band. The pixel size in 
the /3-band images varies in the range between 0.36 and 
0.53 arcsec/pixel, whereas in the //-band images it varies 
between 1.11 and 1.50. In our sample of 152 galaxies, we 
find 102 galaxies with bars in the //-band, and 67 in the 
/3-band. These bars are detected from studying the galaxy's 
Fourier density amplitude profiles. We describe the details 
in the next section. 



3 METHOD AND MEASUREMENTS 

We use a gravitati o nal b ar torque method described by 
[Laurikainen fc Said (|2002D . Bar strengths are estimated 
from the magnitude of the non-axisymmetric gravitational 
perturbation compared to the mean axisymmetric radial 
force field. The relative strength of the perturbation at each 
radius r is defined as the ratio of the maximum tangential 
force Ff^'^^{r) to the azimuthally averaged mean radial force 
FrM: 



(r) 



(1) 



{FR{r)) ■ 

To get a single measure of bar strength, we use the maximum 
of the (3T(r)-profile in the bar region, denoted by Qg. The 
radius where this maximum occurs, is denoted by rg^ . The 
method uses an exponential model for the vertical density 
profile and an empirical relation between the morpholog- 
ical typ e and the radi al to vertical scale parameter ratio, 
hn/zo, (|de Griid IT99i )_ to estimate the vertical thickness. 
The refined method ( Laurikainen et aDl2004bh also allows a 
bulge correction, in which case the different 3D distribution 
of bulge and disk components is taken into account in the 
calculation of the forces. 

To study bar strengths with as little interference from 
the spiral arms as possible, iButa. Block fc Knaoenl (|2003h 
developed a technique for separating the bar and spiral con- 
tributions to the Fourier density amplitude profiles. In this 
study no such bar-spiral separation is made, mainly because 
it is more difficult to apply the method at high redshift. 
Also, the effect of this correction turns out to be small with 
resp ect to the m ean Qg-values in the Hubble type-index bins 
(La urikainen et al..,2007 ) . 

We measure Qg both using the B and //-band images, 
and the results are collected in Table [T] Note that the table 
contains only the barred galaxies from our sample. The clas- 
sification is based on Fourier analysis made in the //-band. 
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Table 1. Barred galaxies in our sample and their bar lengths and bar torque properties. The 
columns are the following: (1) NGC number; (2) Revised morphological type index from RC3; 
(3) r'K2o radius from 2MASS catalogs; (4)-(5) bar lengths estimated by Fourier analy s is in B - 
and H-band, respectively. The latter measurements are taken from L aurikainen et al. 
(6) and (8) bar strength Qg, measured using the H and B- bands; (7) and (9) the radius of Qg in 
the same bands; (10) Qy-value from the B-band, corresponding to rQ^{H). 



Galaxy 


T 




barlen(B) 


barlen{H) 


Qg(H) 




Qa(B) 




Qg(BH) 


(1) 


(2) 


(") (3) 


(") (4) 


(") (5) 


(6) 


{") (7) 


(8) 


{") (9) 


(10) 


JNGO 0150 


3 


56.3 


25 


27 


0.47 


24 


0.55 


25 


0.55 


NGC 0210 


3 


54.1 




46 


0.08 


34 


0.10 


33 


0.10 


NGC 0289 


4 


65.8 




20 


0.21 


13 


0.45 


3.7 


0.30 


NGC 0428 


9 


42.5 


44 


45 


0.24 


44 


0.39 


44 


0.39 


NGC 0578 


5 




18 


20 


0.21 


8.1 


0.26 


7.8 


0.26 


NGC 0613 


4 


118.2 


98 


104 


0.40 


68 


0.53 


72 


0.52 


NGC 0685 


5 


66.7 


21 


21 


0.46 


3.5 


0.49 


9.4 


0.43 


JNGC 0864 


5 


61.8 




26 


0.36 


20 


0.48 


21 


0.48 


NGC 1073 


5 


57.7 


37 


38 


0.62 


23 


0.75 


19 


0.68 


NGC 1087 


5 


59.5 


10 


18 


0.59 


4.5 


0.76 


2.4 


0.68 


NGC 1187 


5 


86.1 




29 


0.21 


36 


0.33 


37 


0.33 


NGC 1241 


3 


56.5 


20 


30 


0.26 


20 


0.30 


24 


0.28 


NGC 1300 


4 


131.8 


85 


87 


0.55 


66 


0.82 


66 


0.82 


NGC 1302 





80.0 


21 


26 


0.08 


22 


0.08 


22 


0.08 


NGC 1317 


1 


63.1 


49 


58 


0.09 


38 


0.10 


35 


0.10 


NGC 1350 


2 


100.1 




81 


0.26 


64 


0.29 


75 


0.28 


NGC 1385 


6 


56.1 


10 


9.3 


0.57 


3.5 


0.64 


5.3 


0.62 


NGC 1493 


6 


66.6 




23 


0.36 


10 


0.63 


11 


0.63 


NGC 1559 


6 


89.7 




17 


0.33 


5.8 


0.52 


10 


0.42 


NGC 1617 


1 


115.7 


16 


22 


0.29 


3.3 


0.35 


1.2 


0.32 


NGC 1637 


5 


76.6 


22 


23 


0.20 


17 


0.18 


15 


0.17 


NGC 1703 


3 


51.4 


10 


11 


0.86 


7.8 


0.10 


10 


0.09 


NGC 1808 


1 


131.1 


86 


87 


0.29 


75 


0.59 


67 


0.55 


NGC 1832 


4 


51.6 


" 


17 


0.20 


12 


0.21 


11 


0.20 


NGC 2139 


6 


43.8 


16 


17 


0.41 


3.3 


0.46 


5.8 


0.38 


NGC 2207 


4 


72.5 


35 


46 


0.33 


29 


0.75 


29 


0.75 


NGC 2442 


4 


126.4 


" 


93 


0.67 


78 


0.99 


76 


0.98 


NGC 2559 


4 


94.7 


27 


33 


0.33 


23 


0.43 


28 


0.41 


NGC 2566 


2 


90.2 


62 


72 


0.33 


52 


0.40 


53 


0.40 


NGC 2964 


4 


58.1 




30 


0.32 


20 


0.32 


20 


0.32 


NGC 3059 


4 


116.6 


20 


20 


0.61 


3.3 


0.59 


5.8 


0.58 


AT/~1<^ 01 /^/^ 

NGC 3166 





81.7 




45 


0.17 


29 


0.15 


25 


0.15 


NGC 3227 


1 


92.6 


82 


75 


0.16 


56 


0.34 


69 


0.30 


AT/~1 00/^-1 

NGC 3261 


3 


61.6 


23 


28 


0.21 


19 


0.27 


21 


0.26 


NGC 3275 


2 


68.7 




41 


0.19 


23 


0.22 


28 


0.19 


NGC 3319 


6 


45.5 


36 


38 


0.55 


14 


0.83 


14 


0.83 


NGC 3338 


5 


91.9 




23 


0.22 


4.5 


0.29 


1.5 


0.20 


NGC 3504 


2 


65.3 


59 


60 


0.30 


29 


0.45 


31 


0.43 


NGC 3507 


3 


88.6 


21 


23 


0.18 


20 


0.20 


19 


0.19 


NGC 3513 


5 


71.7 


27 


28 


0.54 


14 


0.81 


16 


0.80 


NGC 3583 


3 


47.6 


22 


23 


0.25 


17 


0.31 


23 


0.28 


NGC 3593 





101.0 




15 


0.10 


17 


0.10 


23 


0.05 


NGC 3675 


3 


120.7 




30 


0.15 


4.5 


0.20 


1.4 


0.14 


NGC 3681 


4 


38.8 


12 


15 


0.22 


4.5 


0.20 


3.4 


0.20 


NGC 3686 


4 


87.2 


15 


18 


0.28 


11 


0.22 


14 


0.20 


NGC 3726 


5 


103.2 




30 


0.21 


26 


0.22 


28 


0.22 


NGC 3887 


4 


78.5 




41 


0.21 


32 


0.24 


36 


0.24 


NGC 4027 


8 


65.5 


16 


20 


0.62 


3.3 


0.64 


1.9 


0.63 


NGC 4051 


4 


102.6 


41 


45 


0.29 


56 


0.33 


48 


0.34 


NGC 4123 


5 


65.9 




53 


0.43 


38 


0.66 


37 


0.66 


NGC 4136 


5 


53.8 




15 


0.13 


11 


0.13 


10 


0.13 


NGC 4145 


7 


62.8 


15 


20 


0.36 


4.5 


0.28 


1.5 


0.29 


NGC 4151 


2 


88.2 


97 


98 


0.13 


65 


0.20 


66 


0.20 


NGC 4293 





117.1 


71 


68 


0.36 


50 


0.40 


42 


0.39 


NGC 4303 


4 


105.8 




30 


0.27 


41 


0.45 


39 


0.45 



4 T. Speltincx, E. Laurikainen and H. Salo 



Table 1. Continued 



Galaxy 


T 


( ) 


oarlen[r>) 
ill \ 

( ) 


oarLen[H ) 

{ ) 


Qg{H) 


III \ 

{ ) 


Qg{B) 


III \ 

( ) 


Qg{BH) 


NGC 4314 


1 


110.1 


71 


75 


0.46 


53 


0.45 


53 


0.45 


NOG 4394 


3 


97.4 


39 


45 


0.28 


32 


0.27 


33 


0.27 


NGG 4450 


2 


112.0 




30 


0.12 


41 


0.11 


40 


0.11 


NGG 4457 





62.1 


43 


45 


0.09 


32 


0.11 


29 


0.10 


NGG 4487 


6 


67.3 


18 


22 


0.17 


7.8 


0.81 


1.2 


0.24 


NGG 4490 


7 


117.4 




18 


0.33 


7.5 


0.54 


2.6 


0.45 


NGG 4496 


9 




22 


26 


0.36 


7.5 


0.55 


11 


0.51 


NGG 4527 


4 


141.0 




83 


0.18 


77 


0.25 


74 


0.24 


NGG 4548 


3 


105.0 


61 


68 


0.36 


53 


0.36 


52 


0.36 


NGG 4579 


3 


120.0 


46 


45 


0.20 


32 


0.22 


33 


0.22 


NGG 4593 


3 


80.6 




61 


0.32 


43 


0.31 


42 


0.31 


NGG 4618 


9 


73.1 


26 


33 


0.36 


11 


0.64 


12 


0.64 


NGG 4643 





92.5 


62 


68 


0.27 


41 


0.27 


41 


0.27 


NGG 4647 


5 


86.0 


11 


12 


0.24 


11 


0.37 


10 


0.37 


NGG 4651 


5 


69.1 


14 


23 


0.17 


4.5 


0.22 


1.8 


0.16 


NGG 4654 


6 


96.4 


15 


20 


0.14 


11 


0.31 


10 


0.31 


NGG 4665 





110.1 


53 


60 


0.27 


38 


0.25 


34 


0.25 


NGG 4691 





71.5 


45 


45 


0.59 


11 


0.79 


13 


0.73 


NGG 4699 


3 


100.9 




14 


0.08 


10 


0.09 


12 


0.09 


NGG 4781 


7 


69.7 




39 


0.44 


14 


0.60 


15 


0.59 


NGG 4900 


5 


64.5 


14 


18 


0.34 


14 


0.59 


18 


0.49 


NGG 4902 


3 


49.1 




22 


0.28 


17 


0.34 


17 


0.34 


NGG 4930 


3 


58.3 


43 


44 


0.21 


32 


0.27 


35 


0.26 


NGG 4939 


4 


62.7 




17 


0.20 


3.5 


0.22 


1.2 


0.17 


NGG 4995 


3 


57.1 


18 


23 


0.28 


20 


0.27 


25 


0.25 


NGG 5005 


4 


130.8 




45 


0.15 


26 


0.10 


37 


0.06 


NGG 5101 





97.8 


66 


70 


0.22 


43 


0.25 


45 


0.25 


NGG 5334 


5 


51.1 


16 


18 


0.36 


11 


0.39 


12 


0.38 


NGG 5483 


5 


68.5 




13 


0.16 


7.8 


0.25 


8.9 


0.25 


NGG 5643 


5 


99.5 




46 


0.42 


34 


0.48 


35 


0.48 


NGG 5701 





74.3 


47 


50 


0.16 


28 


0.17 


30 


0.17 


NGG 5713 


4 


54.0 




30 


0.43 


4.5 


0.77 


6.0 


0.66 


NGG 5850 


3 


79.2 


74 


90 


0.33 


59 


0.36 


62 


0.36 


NGG 5921 


4 


75.4 


56 


53 


0.42 


47 


0.50 


52 


0.49 


NGG 5962 


5 


49.6 


14 


15 


0.13 


14 


0.12 


12 


0.10 


NGG 6221 


5 




35 


41 


0.44 


29 


0.52 


28 


0.52 


NGG 6300 


3 


116.8 


53 


46 


0.19 


34 


0.29 


32 


0.29 


NGG 6384 


4 


99.6 




33 


0.11 


17 


0.10 


18 


0.10 


NGG 6782 


1 


45.4 


45 


46 


0.17 


24 


0.22 


25 


0.22 


NGG 6902 


3 


53.0 




17 


0.08 


3.5 


0.09 


1.2 


0.06 


NGG 7418 


6 


66.2 




15 


0.20 


15 


0.29 


18 


0.26 


NGG 7479 


5 


87.9 


59 


60 


0.70 


41 


0.88 


37 


0.83 


NGG 7552 


2 


76.4 


74 


70 


0.42 


45 


0.64 


45 


0.64 


NGG 7582 


2 


115.2 


82 


93 


0.44 


55 


0.52 


52 


0.52 


NGG 7723 


3 


68.4 


20 


24 


0.35 


17 


0.28 


17 


0.28 


NGG 7727 


1 


63.2 


15 


27 


0.16 


4.5 


0.16 


3.4 


0.15 


NGG 7741 


6 


55.2 


44 


53 


0.72 


11 


0.96 


15 


0.88 



which classifications are given in iLaurikainen. Salo fc Butal 
l|2004al ). Here no bulge correction was applied, mainly be- 
cause doing such a correction at high redshift is not al- 
ways straightforward. Therefore, the effect of possible ar- 
tificial bulge stretching while deprojecting the images was 
controlled by excluding possible central peaks of Qr-values 
in the bulge dominated region. In the //-band we compared 
our Qg-values withou t applying the bu l ge corr ection with 
the values obtained bv lLaurikainen et al.l (|2004lJ ) where that 
correction was applied. We found that although a bulge cor- 
rection might be important for some individual galaxies, it 



does not affect the mean Qg-va lues in the Hubble type-index 
bins (|Laurikainen et al]|2004bh . 

A similar test was made also in the i3-band, by mea- 
suring the bulge-corrected Qg-values for 26 randomly cho- 
sen galaxies. The bulge correction was ma de by applying the 
2D multi-component decomposition cod e l|Laurikainen et aU 
l2004bl : ILaurikainen. Salo fc Butall2005l '). The method uses a 
Sersic function for the bulge, an exponential disk, and a 
Ferrers or Sersic function for the bar/oval. In order to ob- 
tain robust bulge and disk parameters, a two-component fit 
was made first. Using this fit for the initial parameters we 
then fitted the bar, which is a tricky task in the /3-band, as 
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part of the bar can be obscured by dust and by recent star 
formation. For this reason, we had to occasionally fix the 
bar radius in the fit, based on visual inspection of the im- 
age. Having found an acceptable decomposition, the derived 
bulge model is then subtracted from the original image, the 
image is deprojected to a face-on orientation, after which the 
bulge is added back and the forces are calculated. We com- 
pared the Qg-values before and after bulge correction, and 
found an average difference of 1%. It is therefore safe to say 
that not making the bulge correction will barely influence 
our results. 

The length of a bar can be estimated by looking at the 
phases of the m = 2 and m = 4 components of the Fourier 
density amplitudes: the phase is assumed to be maintained 
nearly constant in the bar region. For the //-band we used 
the bar lengths derived by iLaurikainen et al] (|2004bl ) . For 
the _B-band we derived them in this study, using a similar 
method, and the resulting values (denoted by barlen) can 
be found in Table [1] However, in the B-band the length 
of the bar could not always be estimated by the Fourier 
method, most probably due to the effects of dust that make 
the interpretation of the phase profile more complicated. 
Nevertheless, even for these cases one can still obtain the 
B-band Qrir) profile and its maximum value in the bar 
region, Qg{B). 

In Fig. [T] we show examples of galaxy images and the 
QT-profiles for different morphological types (figures for the 
rest of the galaxies are available in the electronic version). 
The left column shows the _B-band images, whereas the right 
column shows the //-band images in the same scale and ori- 
entation. The center column shows the Qr-profiles in both 
bands (full and dashed line for B and H, respectively) and 
the bar lengths are indicated by a vertical line. It is clear 
from this figure that the _B-band profile is generally less 
smooth due to the influence of the spiral arms. For the early- 
type galaxies, particularly for those with strong bars, the bar 
torques derived from the B and //-band images are gener- 
ally very similar, whereas for later type galaxies the i3-band 
images seem to give slightly stronger bar strengths. 



4 THE BAR TORQUE METHOD AT HIGH 
REDSHIFT 

4.1 Comparing bar strengths 

In order to compare bar strengths in B and //-bands we use 
four approaches, all based on the gravitational bar torque 
method. They are graphically illustrated in Fig. [2] for one 
typical case (NGC 3513). The main idea is to compare bar 
strengths derived from the B and //-band images, either us- 
ing a single number for the bar strength (maximum of Qt{t) 
in the bar region), or by integrating over the area under 
the Qt (T')-profile. The ratio of the obtained bar strengths 
in the two wavelengths is then denoted by Q^y^, where 
i = 1, 2, 3, 4 for the different cases. 

In the first case the averages are calculated using the 
limiting integration radius of rK2(|B) corresponding to the 
radial distance where the /('s-band surface brightness is 20 

^ We obtained the values for »'K2o from the 2 Micron All Sky 
Survey (2MASS) 



mag/arcsec^ . The ratio {Qt{B)) / {Qt{H)) will be denoted 
by Q^y^- In the second approach we use a single number, 
Qg, to measure the bar strengths in the //-band. In the B- 
band we take the Qr-value at the same radius (the ratio 
of Qt's is denoted by Q^gj^)- The third approach is quite 
similar to the second one, but differs in the fact that in 
the i3-band we now take the real Qt maximum closest to 
the radius where the maximum Qt in the //-band appears 
(the ratio is denoted by Qg^/jj)- Our fourth approach uses 
again the area below the Qr-profiles. This time, though, we 
take the integration cut-off at the radius corresponding to 
the length of the bar as estimated from the //-band image 
(the ratio is denoted by Qgjjj)- Compared to Qgjjj, the ad- 
vantage is the elimination of the possible infiuence of spiral 
arms. 

In FigureOwe show how the Q^g/jj's vary with the mor- 
phological type. These plots are based on 99 barred galaxies, 
omitting the three galaxies for which no A's-band luminosi- 
ties were found from 2MASS. For Q^'^there is a slight in- 
crease as we move to later types {Q^g^^jj= 1.1 — 1.7 over types 
through 6, with an average of 1.37), most probably related 
to the prominence of spiral arms well outside the bar region 
in those types. For the other three approaches, we get a more 
constant behavior in the Hubble sequence. Since from type 
7 onward we have only 1 or 2 galaxies per type-index bin, 
statistics there are not reliable. 

Using all galaxies in which a bar was identified in the H- 
band, but limiting to the morphological types T = — 6, we 
find an average bar strength ratio Qg/ff= 1.20, 1.26 and 1.30 
for 1 = 2,3 and 4, respectively. These values are indicated 
with a dashed line in Fig. |3] If we take into account only 
those 64 galaxies for which bars were found in both wave- 

(2) (3) 

lengths, Q^^j:^changes from 1.20 to 1.22, while Q^^^and 

(5^'^(and also Q^gjfj) remain the same. In conclusion, it is 
safe to say that the Qb/h ratio is fairly independent of the 
morphological type, so that we can adopt an average value 
of Qb/h= 1-25 for converting bar strengths between B and 
//-bands. Therefore, while measuring band-shifted optical 
bar strengths in galaxies at z ^ 1 using near-IR images, 
corresponding true near-IR bar strengths can be estimated 
for all morphological types by applying this correction. This 
makes it possible to directly compare bar strengths at low 
and at high redshifts. 

Considering that bars are typically easier to detect in 
the //-band, it may be surprising to find larger bar strengths 
in the i3-band. In principle, this might indicate that the 
density contrast of bars is bigger in the /3-band compared 
to that in the //-band, leading to larger apparent tangential 
forces. However, comparison of the Fourier density ampli- 
tudes indicates that the density variations associated with 
the bars are about the same {{A2 /A2) ~ 1.02). 

Another possibility is that the bulge dilution effect is 
smaller in the i3-band compared to the //-band, so that the 
apparent radial forces are reduced. To investigate this fur- 
ther, we make a crude analytical estimate of the bar strength 
ratio. We approximate the tangential force as Ft oc A2d, 
where A2 is the Fourier density amplitude of the m = 2 
term and d is the disk mass. The radial force is proportional 
to the bulge and the disk mass, but both don't necessarily 
contribute a similar amount: 
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Figure 1. A journey through the morphological types. The left column shows the B-band images, while the right column shows the 
i/-band images. The center column shows the Qy-profiles for the galaxies for B and i/-band in full and dashed lines respectively. The 
bar lengths in both bands are marked with a vertical line. 
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Figure 2. Visualization of the four methods for NGC 3513. From top to bottom: Qg^^^i where we use the area below the Qt curves up 



to tk 



-,(2) 



where we use the maximum value of Qy in the bar region, i.e. Qg, of the i?-band image and the value of the B-band 

-,(3) 



image at the same radius; Q g^^, similar to the second method, but wc now take the maximum value in the B-band image in the 

bar region as well; Qg'^i where we use the area below the curves, taking the bar length in the H-band image to be the maximum 
radius. 



(2) 



where is a radial force, with subscripts b and d for bulge 
and disk, respectively. The mass of the bulge component 
is denoted by b. Here, / is a factor of proportionality that 
depends on the distance from the galactic center (i.e. the 
distance where the Qt maximum occurs). This way, we can 
approximate the bar strength ratio, calculated in B and H- 
bands by: 



iB/H 



Ql 

Qt 



d" + fb^ 



dB + /fei 



dHA!^ 



(3) 
(4) 



(5) 



In case the bulge-to-disk mass ratio &/ci would be the same in 
both bands, the second factor on the right in Eq. ((5} becomes 
unity, and the bar strength ratio simply equals the ratio of 
Fourier density amplitudes. We show the relation between 
Qb/h and A2 /A2 in the upper panel of Fig. [l] It is clear 
from this plot that we cannot simply ignore the bulges. Since 
we do not have (b/d)^ for all our sample galaxies, we make 
a further approximation like: 



Ql 

Qt 



1 + ./ 



{b/d)t 



il)' 



(6) 



where the brackets denote a median, taken from the 26 
galaxies for which we have decompositions in both bands. 
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Figure 3. Q^^^^versus morphological type. For Q^j^we see a slight increase as we go to later types, which is an effect of the prominence 

of spiral arms in those types. For the other three methods, the Qg^jj's are fairly independent of the type except for T > 6, where the 
number of galaxies and hence the reliability of the statistics drop considerably. The dashed line marks the average value through those 
seven types. The length of the vertical bars indicates the standard deviation for each type bin. 



For this median, we find a value of ~ 0.41, consistent with 
bulges being red with respect to disks. To calculate an ap- 
propriate value for /, we approximate the bulge by a point- 
mass, while the disk is represented by an exponential with 
scale length Hr. Note that rg^ occurs on average at about 
one exponential scale length: for r/hn « 1 we have / ~ 3. 
The relation from Eq. ((6)l with this value for / is shown in 
the bottom panel of Fig. |4] There is a clear improvement, 
confirming the fact that Qb/h> 1 is mostly due to reduced 
bulge dilution in the _B-band. 

The bar lengths in the two wavelengths are very similar, 
the one in the B-band usually being the slightly shorter one. 
The left panel of Fig.[5]shows a tight correlation between the 
bar lengths in the two bands, the lengths being normalized to 
r_R-2o- The radius of the maximum bar Qt, fQg, is well corre- 
lated with the length of the bar (estimated from the Fourier 
phases), in agreement with 'Laurikain en. Salo fc RautiainenI 
In the right panel of Fig. [5] we show rq^ relative 
to the bar length, compared in the two bands. The outlying 
points are due to underestimated bar lengths in the B-band, 
most likely because part of the bar is obscured by dust. On 
the other hand, the maximum of the Qr-profile still appears 



at the 'right' place, causing the ratio of these two to be over- 
estimated. 

4.2 Estimating za through rx2o 

The gravitational bar torque method makes explicit use of 
the radial scale length of the disk, hn, to obtain th e vertical 
scale height, zo, through the empirical relation by Ide GriisI 
(| 19981 ') ■ where zq/Hr depends of the morphological type. We 
investigate here whether we can also get similar results with- 
out using this relation, by estimating the vertical scale height 
by means of the rK2o radius. We calculate the bar strengths 
again for the galaxies in our sample using the following ap- 
proximation for zq: 

zo / ^0 \ 
zo = rKso ~ ( ) rK2„ (7) 

Here the ratio {zo/rK^o) is calculat ed using the original 
Zo values estimated via the Ide GriisI |l998|) relation. With 
this approximation, we make thin disks (mainly in late-type 
galaxies) a bit thicker, and thick disks (mainly in early-type 
galaxies) a bit thinner. As an average value for all morpho- 
logical types in our sample we find {zo/rif^a) ~ 0.08. This is 
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Figure 4. The influence of bulge dilution. The top panel shows the relation between the bar strength ratio and the density contrast. 
The dashed lines indicate the median values. The bottom panel takes into account the different bulge-to-disk ratios in the two bands. 
For a full description, see Section l4.1l 
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Figure 5. The left panel of this figure shows a tight correlation between the bar lengths, scaled by tk2qi in both bands. The right panel 
shows the correlation between tq^ compared to the bar length, in both bands. The dashed lines indicate the unit slope. 



by no means a universal value, but could be used for spirals. 
If one was to study SO galaxies, another value would most 
likely be more suitable. In Fig.[6]we show how much the bar 
strength is affected when using this approximation. In both 
bands we compare the ratio of bar strengths derived using 
the above approximation, and the original method where 
the radial scale length was used for estimating the vertical 
thickness. There is a clear type dependence visible, which is 
no surprise since zo/Hr depends on the morphological type. 
Thus, estimating through rK2o implies that bar strengths 
for the early type galaxies would be overestimated, whereas 
for late type galaxies they would be underestimated, but 



by no more than 10-15%. So, within these errors, the Qg- 
method is still useful even if the morphological type and/or 
radial scale length can not be determined. Also note that 
the effect of this approximation is similar in both bands. 



4.3 The bar torque method using degraded 
images 

So far we have tested our method on images from a local 
galaxy sample. In order to test whether the method is valid 
also for images with lower resolution, like those obtained at 
higher redshifts, we degraded the images of our sample. As 
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Figure 6. Change of Qg values when zq is estimated from Tfc^^ rather than /i/j. A clear type dependence is noticeable. For early types 
the bar strengths will be overestimated, whereas for late types they will be underestimated. 



this is not supposed to be a test for any particular telescope 
or survey, we don't follow the d etailed technique as used e.g. 
by Ivan den Bergh et al.l (|2002l . and references therein). We 
use the nearest neighbour algorithm, i.e. we take 2" pixels 
of the original image, and put the information from one of 
those, the nearest neighbour, to the resulting pixel in the 
degraded image. This way we keep the same noise level in 
the degraded images. There is no correction made for seeing. 
Since for near-IR observations of high redshift galaxies the 
rest-frame wavelength is in the optical, the image degrada- 
tion was made for the _B-band images. 

In Fig. [7] we show the change in bar strength Qg after 
degrading the images by different factors. After a degrada- 
tion by a factor of 2 (top left panel) or 4 (top right panel), 
the Qg values barely differ from the original ones, indicated 
by the tight correlation in Qg obtained from the original 
and the degraded images. Using a factor of 8 the correla- 
tion still holds (bottom left panel), but using a factor of 16 
only a small percentage of the galaxies with very low bar 
strengths still give nearly the same Qg as before the degra- 
dation (bottom right panel). In Fig.|8]we show an example of 
an extreme case, NGC 6221, for which the whole Qr-profile 
stayed almost unchanged even after degrading the image by 
a factor of 16. 

In Fig. [9] we translate all this into pixel sizes. For each 
barred galaxy in the B-band we calculated the bar length in 
pixels and the deviation from the original bar strength after 
the degradations by a factor of 2, 4, 8 and 16. Therefore, 
each of those galaxies has four corresponding points in the 
plot. In general, as long as the total length of the bar is 
at least 10 pixels, the difference between the bar strength 
derived from the original and the degraded image remains 
typically less than 10% (indicated by the dashed lines in 
the figue). Thus, the Qg method can estimate bar strengths 
quite well, even from low resolution images. With the 0'.'05 
pixels from the ACS, this means that at z = 1, for any bar 
that is at least 2 — 3 kpc (depending on the cosmological 
model), we can reliably estimate its strength. 



5 DISCUSSION 

Estimating bar strengths, even in the local Universe, has 
proved to be far from trivial. The main parameters used 
to characterize bar strength are the ellipticity of the bar e, 
the bar torque Qg, as well as the length and the relative 
luminosity of the bar. These parameters do not show sim- 
ilar trends in the Hubble sequence, mainly because e, and 
particularly Qg are diluted by the axisymmetric force field, 
largely re lated to the more massiv e bulge s in the early-type 
galaxies (|Laurikainen et al]|2004al . |2007| ). This means that 
massive bars (A2 amplitude is large, and/or long compared 
to the scale length) do not necessarily have the strongest bar 
torques. For the characterization of bars at different redshifts 
all these parameters are expected to be important. This is 
also one of the reasons that makes our bar strength calcula- 
tion particularly useful: since our method utilizes azimuthal 
Fourier decomposition of intensity, we obtain, besides Qg, 
also an estimate for the relative luminosity of the bar in 
terms of Fourier density amplitudes. Moreover, the length 
of bar can be estimated based on the phase of the m — 2 
and m — 4 Fourier amplitudes. 

Questions have been raised, though, whether calcu- 
lating bar strengths is feasible at high redshifts . Namely, 
it has been argued by iMarinova fc Joged (|2007l ) that the 
bar torque method probably fails because (i) it is not 
suitable for noisy images, (ii) it is not applicable for im- 
ages with low resolution, and (iii) it makes explicit use 
of the vertical scale height and morphological type, which 
can not always be reliably derived for high redshift galax- 
ies. However, it is important to note that two types of 
bar strength methods have been used in the lite rature: 
the Cartesi an method (|C|)uillen^ Frogel fc Gonzalez! 



1994 ). 



2003), 



applied by iButa fc Blo ck ('2001) and 'Block et al.| 
and the polar m ethod dLaurik aincn fc Salo 200 2: see als o 
Salo et aTl 19991) . applied by iLaurikainen et a.l] (l2004a ). 
Laurikainen et all (|2004bl . bOOel ) and lButa et all (|2004l2005l . 
2006), and also used in the present study. The Cartesian 
method is indeed very sensitive to pixel-to-pixel noise, easily 
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Figure 7. Correlation between the Qg values after degradation by different factors and the original values. Up to a factor of 8 the 
correlation is very good. 



leading to spurious peaks in the force ratio, particularly in 
the outer parts of the images where the S/N-ratio is small. 
However, the polar method is not very sensitive to noise, 
because the use of a Fourier density decomposition prior 
to the calculation of the forces implies averaging over the 
azimuthal angle (for a comparison between the p o lar an d 
Cartesian methods, see ISalo. Laurikainen fc Butal l|2004 )). 
In fact, since the calculation of bar torques integrates over 
the whole luminosity distribution, the force evaluation (with 
the polar method) can be expected to be even less sensitive 
to noise than the method where ellipses are fitted to the in- 
tensity contours, averaged over very narrow radial bins. The 
general characteris tics of the Cartesian and p olar methods 
were compared by [Laurikainen fc Said (I2OO2I ): the robust- 
ness of the polar method with respe ct to noise is also strik- 
ingly illustrated in lSalo et all l|2004 ). 

The two other issues mentioned by iMarinova fc Joge3 
have been address e d in this study, and partly also by 
[Laurikainen fc Said (|2002f ). We found that our Qg-method 
is not overly sensitive to image resolution; the local OSUBSGS 
galaxies observed with the pixel size of 0'.'44 can be degraded 



by a factor of 8, and we still obtain an indistinguishable av- 
erage bar strength. It seems that the method can be applied 
within a 10% accuracy, as long as the total bar length is at 
least about 8 pixels. A s a measure of comparison, we note 
that Ijogee et al.l (|2004l ) find bars in their gems survey with 
a semi-major axis length of 0'.'15 - 2"2. The ACS has 0'.'05 
pixels, so these bars have total lengths ranging from 6 — 88 
pixels. Therefore we conclude that the Qg-method would 
give a reliable estimate of the bar strength for most of those 
galaxies. As stated above, for a bar length of 8 pixels we find 
an average error of 10%, and this decreases to a mere 1% 
when the bar length is at least 40 pixels. 

The estimation of the vertical scale height in the bar 
torque method might appear problematic if the radial scale 
length of the disk is not well known. Probably no detailed 
2D multi-component decompositions can be applied at high 
redshifts. However, a one dimensional bulge/disk decom- 
position using an azimuthally averaged image is even bet- 
ter to estimate Hr i n a re liable manner. It was shown 
by [Laurikainen et al] l|2006l ) that the ID method gives 
very similar /i_H-values to the more sophisticated 2D multi- 



Comparison of bar strengths in optical and near-infrared 13 



1000 



800 



600 



400 



200 



original image 



B-band, NGC6221 




200 400 600 800 1000 



500 
400 
300 
200 

100 





degrade 2 



4 



100 200 300 400 500 



252 
126 




degrade 4 



62 124 




a 0.2- 



126 252 



400 



r (pixels) 



Figure 8. Resulting Qt profiles for NGC 6221 after degradation are shown on the right with the radius in pixels according to the 
original image. On the left we show the corresponding images. 



component approach. In order to estimate the vertical scale 
height, we should also be able to classify the morphological 
type of th e galaxy. Howev er, in order to use the empirical 
relation of Ide Gr"ii3 (|l998l ). it is enough to divide the disk 
galaxies into 3 different Hubble type bins. Also, no ad hoc 
knowledge of the barred properties is required. In case we 
have no information of the morphological type of the galaxy, 
we can use the ^0/^X20 ratio to estimate the vertical thick- 
ness, in which case we can still measu re the bar strength 
within ± 15%. Finally, it was shown by iLaurikainen fc Said 
that the bar torque method is practically indepen- 
dent of the model used for the vertical distribution of the 
light, as long as the scale height is fixed; therefore for sim- 
plicity an exponential function has generally been used. 



Based on the above discussion we believe that the 
bar torque approach is promising for characterizing bars 
even at high redshifts, providing useful complementary in- 
formation to the bar ellipticity. Both the ellipticity and 
the bar torque are affected by the underlying axisymmet- 
ric potential, usually the bulge. However, the bar torque 
seems to be a better indicator of the systematic dif- 
ferences between bars in the Hubble sequence: Qq in - 
creases toward later types l|Laurikainen et al]|2004bl , l2007h . 
whereas e is maintained nearly constant in all Hubble types 



jMenendez-Delmestre et al] l2007l : ILaurikainen et al" I l2007l : 
i Marinova fc Jogoc 2007). 

The Fourier method is also very useful for identify- 
ing bars in galaxi e s, as sh own for the OSUBSGS galaxies by 
ILaurikainen et al] (|2004ah . In that study several different 
types of diagrams were used for the identification, but the 
most important criteria were to detect the m = 2 Fourier 
amplitude above the noise in the bar region, and to look at 
the phases of the m = 2 and m = 4 Fourier components: 
the phases are maintained nearly constant in the bar re- 
gion. This approach is expected to pick up classical bars, 
but not the more oval-like or spiral-like bars, generally clas- 
sified as SAB types in RC3. The second criterion is actually 
the same as in the approach where ellipticities are used to 
identify bars: besides requiring that the ellipticity increases 
steadily to a global maximum, the position angle is also as- 
sumed to remain constant in the bar region. Therefore it 
is not unexpected that similar bar frequencies in the OS- 
UBSGS sam ple have been obt ained by applying the Fourier 
method by Laurikainen ot al. (2004a"), and by the ellipticity 
approach by Marinova fc Jogcc (2007). Using an inclination 
limit of 60°, the found bar frequencies are 62% and 60%, 
respectively. The ellipticity approach has been used also by 
iMenendez-Delmestre et ail (|2007l ) , applied for a magnitude- 
limited sample of 2MASS galaxies: using an inclination limit 
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Figure 9. Error in per cent in relation to tlie total bar length in pixels. Individual galaxies are plotted with squares. Each galaxy is 
represented by four data points, corresponding to different degradation factors. The full line shows the average error, using logarithmic 
bins for the pixel size. It is clear that for bars that occupy 10 pixels or more, the average error will be below 10% (indicated by the 
dashed lines). 



of 65°they found that 59% of the galaxies are barred, which 
is exactly the same value as found by us using the Fourier 
approa ch with the same inclination limit (Laurikai nen et al.l 
l2004bl ). In all these studies the number of bars increases 
by some 10% if the criterion of the Fourier phases/position 
angles is omitted. It seems that both methods pick up the 
classical bars well, with results slightly larger t han the 56% 
of SB-type bars visually identified by Eskridg e et alj (|2002l ) 
in the near-IR for the OSUBSGS sample. This is good news 
because it means that both Fourier and ellipticity methods 
can be used to estimate bar frequencies at high redshifts. 
However, it is worth noticing that the above comparisons 
have been made for spiral galaxies, which generally have 
small bulges. If the bulges are large, which is the case for 
some early-type disk galaxies, the bars might be overshad- 
owed by the light of the bulge in the surface brightness 
profile. In that case, due to the azimuthal averaging, the 
Fourier method is expected to be more sensitive to pick up 
the non-axisymmetric component, even if that structure is 
not prominent in the contours of the two-dimensional surface 
brightness distribution, used to identify bars in the elliptic- 
ity approach. 

Probably the most unexpected result from this paper is 
that the bar strengths in the B-band are larger than in the 
ii"-band, especially since bars are found to be more promi- 
nent in the infrared (e.g. lEskridge et al.l [20021 ) . We found 
(Section 14. ip that different bulge dilution in optical and in- 
frared explains the larger bar strengths in the optical. The 



bulge is less prominent in the optical, underestimating the 
radial force due to the bulge and hence artificially enhancing 
the bar strength. 



6 SUMMARY AND CONCLUSIONS 

We have used a gravitational bar torque method to compare 
bar strengths derived from the B and ff-band images for the 
galaxies in the OSUBSGS sample. We presented four different 
ways to do the comparison. Further we subjected the method 
to various tests simulating a high redshift environment. We 
estimate zo through rK2o, and degrade the B-band images 
using a nearest neighbour algorithm before calculating the 
bar strengths. Our main conclusions can be summarized as 
follows; 

(i) Quite unexpectedly, we found that bars appear to 
have stronger torques in the optical than in the near-IR: the 
conversion factor between B and il/-bands is 1.25. We also 
showed that this is due to reduced dilution of radial forces 
by relatively smaller bulges in the optical. 

(ii) The scaled bar lengths, barlen/rK2o, and 
rQg /barlen are similar in the B and Jf-bands, with the bars 
in the B-band being on average slightly shorter than bars in 
the H-hand. 

(iii) An approximation of the vertical scale height while 
calculating Qg can be estimated from rjfjg. Resulting bar 
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strengths will be overestimated for early types and underes- 
timated for late types, by less than 15%. 

(iv) The Qg method is found to work well even on low 
resolution images. We degraded our B-band images using 
the nearest neighbour algorithm and find that as long as 
the total bar length is at least ~ 10 pixels, the resulting Qg 
value is typically within 10% of the original one. 

We have shown that the Qg method provides a reliable 
tool for measuring bar strengths at high redshifts: with the 
pixel resolution of 0'.'05 of the ACS, bars with 2 — 3 kpc 
can still be reliably measured at z ~ 1.0. 
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